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The dynamics of granular materials, mostly radial and axial segregation in horizontal rotating 
cylinders filled with millet and poppy seeds, is studied by Magnetic Resonance Imaging (MRI). For 
the first time, full 3D structures and real-time 2D MRI movies showing the progress of segregation 
over many hours are reported. Data are acquired with sufficiently high temporal and in-plane spatial 
resolutions (74 ms and 0.94 mm x 0.94 mm, respectively), giving new insights into the underlying 
mechanisms. The millet and poppy mixture composition is calibrated based on the signal intensity 
and is quantified throughout segregation. As for radial segregation, millet and poppy mixture core 
formation is observed in cylinders of 75% and 82% filling level. The size of the core is calculated 
and the avalanche layer thickness is therefore determined. 2D MRI movies showing real-time radial 
segregation suggests that initial radial segregation is due to the stopping and percolation of poppy 
seeds. Axial segregation is characterized by the formation, traveling and merging of poppy-rich 
bands. In all cases studied, the formation of poppy-rich bands is observed, after which individual 
bands start to travel at ^3 /im s _1 until they are within ^3 cm of a stationary band. Adjacent 
bands then merge into a single, enlarged poppy band as millet seeds move out of the merging region. 
In both radial and axial segregation, core diffusion is shown to be another mechanism pathway for 
segregation besides the free surface. 

PACS numbers: 45.70.Mg, 76.60.Pc, 45.70.n, 45.70.Ht 



I. INTRODUCTION 

Granular materials, including grains, gravels, and 
many powders, are aggregates of large solid particles. 
They are known to have a highly complex disordered 
structure, nonlinear internal friction, tend to form stable 
arches and voids, and inelastic collision between parti- 
cles under load [1]. As a consequence, granular materials 
possess many physical properties which are not observed 
in classical solid or liquid state physics, such as internal 
stress fluctuations [2], intermittent avalanches [3], segre- 
gation [4], and spatial pattern creation j5(]. 

Horizontal rotating drums are common in industry, 
such as tumbling mixers and reactors, and the segrega- 
tion of granular materials in this geometry is of much 
interest. In such systems, segregation can occur radially 
and axially. This is a well-known phenomenon and has 
been reported many times in the literature. However, 
there are still areas of controversies, such as the mecha- 
nism of radial segregation, the role of core diffusion ver- 
sus that of the free surface, which will be reviewed in 
turn below. In addition, this work also reports new ob- 
servations of radial segregation and band traveling and 
merging. This area has only been touched upon in the 
literature, so a review of the current knowledge is also 
included. 

The mechanism that drives radial segregation has been 
considered in the literature. There is an argument that 
radial segregation is due to competing attractors, one 
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near the outer edge that has higher affinity for larger 
particles and one near the rotation axis that has higher 
affinity for smaller particles |6]. However, a different ap- 
proach is based on the statistical mechanics in which the 
granular material is treated as a sample of many small 
particles and their collective behavior depends on the 
sample volume and compactivity [7]. Another explana- 
tion attributes this to the percolation and stopping of 
small particles [8]. It is reasoned that as a mixture of 
particles flow down the slope in the avalanche layer, the 
smaller particles carry less momentum due to their lower 
mass and therefore are more likely to stop in the middle 
of the slope during collisions between particles while the 
larger still keep rolling downward and get into the outer 
layer. This creates segregated layers of smaller particles 
near the bottom of the flowing layer and larger near the 
top [5j]. As a result, smaller particles do not tend to reach 
the cylinder wall at the bottom of the inclined avalanche 
region but instead form a core along the rotation axis. 
This process is very rapid, which renders radial segrega- 
tion visible within one revolution and possibly complete 
within a few. 

The free surface has mostly been thought to play the 
sole role in both radial and axial segregation so far, espe- 
cially in modeling work. Lai et al. simulated radial seg- 
regation due to friction by using the BTW (Bak, Tang, 
and Wiesenfeld) sand-pile model [9]. However, the ro- 
tating cylinder was assumed to mimic a sand pile, so no 
consideration was given to the geometry of the cylinder, 
nor its rotational speed. These details were, however, in- 
cluded in a later and more detailed study, together with 
the dynamics of the particles in the thin avalanche layer 
[lot . This was based on the density difference of the same 
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sized-particles and the mechanism of the segregation was 
assumed to be purely due to the avalanche layer. As will 
be shown later, this is only valid for systems of less than 
75% filling. Other 2D models have also been considered 
but restricted to either single particle dynamics p, or 
low rotational speed [ll|. In all cases, the avalanche layer 
was considered to be the only mechanism for segregation. 
Other authors propose that diffusion inside the core may 
play a role [12[- From geometrical considerations, it is 
evident that all seeds in the bulk reach, and are there- 
fore influenced by, the avalanche layer in every rotation 
in a cylinder less than 50% full. However, at higher filling 
fraction, there is a cylindrical core about the rotation axis 
which may never reach the avalanche layer. It is in our 
interest to study such systems to examine the underlying 
mechanisms. 

Upon radial segregation, further rotations lead to the 
formation of axial bands at equal spacings, as reported 
by some authors to date [l2[. However, the evolution of 
bands does not stop at band formation but develops into 
band traveling and merging over a longer time scale. In 
one study, the particle velocity in the flowing layer was 
found to be higher for larger particles, and thus thought 
to be the main drive for axial segregation [13] using the 
camera technique. This method involves the imaging of 
the cylinder surface by a charged coupled device (CCD) 
camera and the velocity is calculated by tracking individ- 
ual particles. The data were then, together with the fric- 
tional properties, successfully implemented into a Monte 
Carlo simulation to reproduce band merging. However, 
this model is mainly qualitative, segregation is modelled 
to occur much faster than reality and has no prediction of 
what will happen beyond the first merging. Using Mag- 
netic Resonance Imaging and camera techniques, other 
authors also reported band merging [HEl. The MRI 
data give some qualitative insight into the progress of 
the merging bands but only in terms of visual observa- 
tion and is not sufficient as only 2D stationary images 
were acquired. A continuum model based on two local 
variables, concentration difference and the dynamic angle 
of repose, was subsequently proposed to fit the data ob- 
tained by the camera technique [15|. However, the model 
was limited to predicting when the merging would occur, 
without revealing what would happen to the particles at 
the front of and between the merging bands. Nonetheless, 
the author was successful at working out a logarithmic 
rate of reduction of the number of bands over time, which 
is again confirmed in a later study [lq |. No real 3D data 
are available to give full description of the mechanisms 
of band merging. 

In this study, radial and axial segregation are investi- 
gated using MRI, which is very useful for the study of 
granular flows as most of the systems are opaque. Dyed 
beads have been used but only surface observations can 
be made and are inappropriate for examining the state of 
segregation inside the bulk material. This problem can 
be rectified by using particles with different MRI activ- 
ities, which renders them distinguishable when mapped 
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FIG. 1: 3D RARE pulse sequence diagram. Homospoil can 
be used before and after the 90° pulse. 
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FIG. 2: 2D FLASH pulse sequence diagram. The readout 
gradient is often prolonged to act as homospoil. 



in space and time. The theory of MRI can be found in 
an excellent text by Ref. [17| . 

The main pulse sequence used in this work is 3D RARE 
(Rapid Acquisition with Relaxation Enhancement) with 
the pulse diagram shown in Fig. [T] [18]. This is a 
fast imaging technique that samples multiple-line k-space 
from each excitation, thus reducing the acquisition time. 

The disadvantages of the RARE sequence, though hav- 
ing good signal to noise ratio (SNR), are that it cannot 
image moving objects and its acquisition time is of the 
order of minutes for 3D or many seconds for 2D images, 
therefore not sufficiently rapid to acquire images in real 
time. For this reason, 2D FLASH (Fast Low Angle SHot) 
is instead employed, the pulse sequence of which is shown 
in Fig. [2] [19]. The main feature of FLASH is the very 
low tip angle used, typically 5° — 10°, and thus has an 
inherently low SNR 

In this work, full 3D structures and real-time 2D 
movies are presented, providing new insights into segre- 
gation phenomena. New observations of radial and axial 
segregation are made and issues brought up in the liter- 
ature review are addressed, by conducting a number of 
experiments listed below: 

a. New observations: core formation throughout ra- 
dial segregation in cylinders of 75% and 82% full 
cylinders. The avalanche layer thickness can then 
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be inferred from the core size. Cylinder length is 
15 cm and 12 cm, respectively, but has no effect on 
radial segregation. 

b. Mechanism of radial segregation: real-time 2D 
movies showing the stopping and percolation mech- 
anism. Experiments are conducted with an 82% 
cylinder, rotated at 10 rpm. Radial segregation oc- 
curs slowly in such a full cylinder at a low rotational 
speed, therefore assisting data acquisition. 

c. The role of core diffusion versus that of the free 
surface: can be seen in both radial and axial segre- 
gation and is examined in both sections IIII Al and 
IIIIB1 Supporting data are acquired during radial 
segregation in cylinders of 75% and 82% full; as well 
as during axial segregation in an 82% full cylinder. 
Only cylinders of more than 75% filling are cho- 
sen because the effect of the free surface movement 
and core diffusion can be uncoupled in a core buried 
deep under the avalanche layer. 

d. New observations: band traveling and merging. 
The long time scale 3D RARE data are acquired 
for a 50% full and 50 cm long cylinder, rotated at 
35 rpm. The short time scale 2D FLASH data are 
acquired for a 50% full cylinder, also rotated at 35 
rpm, with the lenth of 12 cm to fit within the FOV. 
Only 50% full cylinder is examined because axial 
dynamics occurs much more rapidly compared to 
that in 75% or 82% full cylinders. 



II. METHODOLOGY 

The system to be studied is a closed horizontal cylinder 
of 48 mm I.D., filled with a mixture of poppy and millet 
seeds. Both seeds are MRI active due to their oil content 
but poppy seeds give a signal (5000 au pixel -1 ) 10 times 
stronger than millet (500 au pixel -1 ). Poppy seeds are 
generally ellipsoidal with a major axis of 1 mm and minor 
axes of 0.7 mm, approximately; whereas millet seeds are 
more or less spherical of 2.5 mm in diameter. Poppy 
and millet have similar density (1100 kg m -3 ) so their 
segregation is expected to be mainly size-driven (2l|. It 
was found from previous experiments that segregation 
was most obvious when the poppy seeds were less than 
millet in volume so a composition of 25%: 75% vol (poppy 
: millet) was adopted for all experiments. The angle of 
repose of poppy seeds is approximately 2° higher than 
that of millet, which is just over 30°. 

The cylinder length is adjustable and chosen depending 
on the application. A typical length is 12 cm so as to fit 
in the FOV of the magnet. However, a 50 cm length 
has been used so that band traveling and merging can 
be observed. The cylinder filling fraction is also varied 
depending on the purpose of the experiment. 75% and 
80% full cylinders are used for experiments studying the 
role of core diffusion versus that of the free surface as 



Magnet 



Motor 



Magnet 



FIG. 3: Schematic diagram of the experimental setup and its 
coordinates. 



these two effects can be uncoupled. 50% is chosen to 
study axial segregation as it is conventional and occurs 
much more rapidly to assist data acquisition. 

The seeds are premixed externally before loaded into 
the Perspex tube, which is rotated by an electric mo- 
tor located about 2 m away from the horizontal Bruker 
Biospec AV magnet (85 MHz, ^-resonance, 2 T) as 
shown in Fig.[3j The two adjustable dividing plugs allow 
variable tube lengths of interest. 

All 3D MRI data are acquired using a RARE sequence 
while the motor is at rest. The RARE pulse sequence is 
implemented as follows: The 90° and 180° soft pulses are 
both 512 ms sine pulses. The spin-lattice relaxation time 
T\ is measured using an inversion recovery pulse sequence 
in which a 180° hard pulse is first applied, followed by 
a variable delay time before a 90° pulse and the signal 
is acquired. T\ is then calculated from the plot of signal 
against delay time. Both poppy and millet seeds have 
rather short 7\, i.e., 200 ms, and only 32, i.e., RARE 
factor, echoes are collected per excitation. The field of 
view (FOV) is 6 cm in the x, 12 cm in y and 3 or 6 cm 
(depending on the system being 50% filled or higher) in 
the z-direction. The resolution is thus 0.94 mm pixel -1 , 
giving SNR = 10 after four averages. The echo time is 5.1 
ms, with a 300 ms repetition time and the whole acqui- 
sition time of just under 4 min for a 128 x 64 x 64 image. 
It is also noted that the magnetisation is returned to 
the origin of k-space after each RARE readout sequence. 
This means that all k-space lines are acquired in the same 
direction and the odd-even echo problem is thus reduced 

Radial and axial segregation are also captured in real 
time by using 2D FLASH. The pulse used is a 256 ms 
Gaussian pulse with an excitation pulse angle of 10°. The 
echo time TE = 1.3 ms and repetition time TR = 4.7 ms. 
As for radial segregation, an 82% full cylinder, rotated at 
10 rpm is chosen. Radial segregation occurs quite rapidly 
so only one average is used to gain temporal resolution. 
In order to gain enough SNR, a 5 cm thick slice in the xz 
plane is imaged due to negligible axial segregation within 
the first revolutions, with the FOV of 6 cm x 6 cm so as 
to avoid fold over. The in-plane resolution is chosen to be 
64 x 64, resulting in temporal resolution of 300 ms. As for 
axial segregation, the fuller the cylinder, the slower seg- 
regation. Therefore, only a 50% full cylinder is imaged so 
as to keep segregation within practically short times, i.e., 
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~ 13.5 min instead of many hours. A 2 mm thick slice in 
the yz plane, running along the rotation axis and cutting 
through the center of the sample, with the FOV of 12 cm 
x 3 cm, resolution of 128 x 32 is imaged, repetition time 
of 150 ms. Due to the thin slice used and the inherent low 
SNR of FLASH, four averages are used to collect enough 
signal in this case. FLASH uses a single gradient echo 
to acquire data, which reduces the repetition time to a 
few milliseconds per k-space line. The main disadvan- 
tage of FLASH is the inherent low SNR since the signal 
is only proportional to sin#. However, the longitudinal 
magnetisation is then proportional to cos and is mostly 
compensated by spin lattice (Ti) relaxation, especially in 
short T\ systems, thus becomes saturated quickly and an 
infinite number of images could, in principle, be acquired 
without losing SNR [S3. 

In all MRI images, the signal intensity in each pixel is 
obtained in terms of an arbitrary unit (au pixel -1 ) com- 
ing from both the poppy and millet seeds and is thus 
proportional to the amount of seeds present. A calibra- 
tion of the signal intensity as a function of the poppy 
concentration by volume is shown in Fig. [4j It can be 
seen that the signal intensity is linearly proportional to 
the poppy concentration except when this falls below 20% 
by volume. This could be attributed to the low signal in- 
tensity of the millet seeds (i.e., 500 au pixel -1 ) compared 
to that of noise (i.e., 195 au pixel -1 ). Nevertheless, the 
local mixture composition in each pixel can be accurately 
quantified directly from the signal intensity. 
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FIG. 4: Calibration of the signal intensity as a function of the 
poppy concentration by volume for the 3D RARE sequence 
proposed. 



All MRI data are processed using an in-house soft- 
ware package. The raw data (after Fourier transform) are 
gated at 195 au pixel -1 to suppress the background noise. 
Segregation is found to be rather symmetrical about the 
rotation axis, and radial segregation in particular is also 
invariant along the cylinder length. Therefore, projected 
images can be used to study radial and axial segregation 
for representative data sets. In addition, signals towards 
the ends of the receiving coil (in the ^/-direction) start 



to drop due to the magnetic field inhomogeneity. There- 
fore probe correction is required for all data sets by using 
correction factors worked out from a phantom filled with 
seeds. No correction is necessary in the x or z direction. 



III. RESULTS AND DISCUSSIONS 
A. Radial segregation 

In this section, we report core formation during radial 
segregation in cylinders of 75% and 82% filling. Based 
on the core size, the avalanche layer thickness can be 
inferred, which is then compared with measured values 
stated in the literature. In addition, the role of core diffu- 
sion is confirmed by examining the poppy concentration 
in the core during radial segregation. Also, the mecha- 
nism of initial radial segregation is addressed by looking 
at real-time 2D FLASH movies. Such full systems have 
been chosen partly due to the slower rate of radial segre- 
gation and partly because core formation occurs. 

Figure [5] shows the 3D structure of a 75% full and 15 
cm long cylinder, rotated at 35 rpm, prior to and after 
960 s (560 revolutions) of rotation. The concentration 
of the poppy seeds is linearly proportional to the signal 
intensity, hence can be easily calculated as shown in sec- 
tion HI The color bar shows the corresponding poppy 
concentration for each color, however for easy identifica- 
tion, red indicates poppy-rich region, green is an equal 
mixture of poppy and millet, cyan is millet, and blue is 
void. It can be seen that seeds are initially well-mixed 
with mostly green color in the cut-through view. Radial 
segregation only comes to completion after 960 s when 
most of the poppy seeds concentrate in a cylindrical core 
along the entire length of the rotation axis, whereas the 
rest of the drum is millet, i.e., cyan color. 

A series of these stationary 3D images are acquired at 
many intermediate rotations showing the progress of ra- 
dial segregation. Each 3D image is then collapsed into 
a 2D projection in the xz plane due to negligible axial 
segregation such as those shown in Fig. [6j Since radial 
segregation concerns mainly the traveling of seeds radi- 
ally and the 2D projection is rather symmetrical about 
its center point, it is possible to comprehensively repre- 
sent the status of radial segregation by a ID profile across 
the cylinder radius. This is done by averaging the poppy 
concentration at each distance from the core center in the 
2D projection to get a ID profile as shown in Fig. [7(a). 
A similar procedure is also done with an 82% full, 12 
cm long cylinder, and results are shown in Fig. [7(b). A 
zero distance to core center indicates the center point of 
rotation, i.e., on the rotation axis; whereas a maximum 
distance, i.e., 30 mm away, means the cylinder wall. 

At both filling fractions, it is found that subject to ro- 
tation, poppy seeds immediately move inward, forming a 
core so that only millet seeds are visible from the outside. 
This can easily be seen in Fig. [6(a) after 2 s of rotation, 
or in Fig. where there is, as soon as rotation starts, an 
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FIG. 5: 3D structure of radial segregation in a 75% full cylinder after (a) s and (b) 960 s of rotations at 35 rpm. Color codes: 
red (poppy), green (mixture of poppy and millet), cyan (millet), and blue (void). FOV x x y x z: 6 cm x 12 cm x 6 cm. The 
coordinates are in mm. 
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FIG. 6: 2D xz projection, produced from the 3D structures in Fig. [5] showing radial segregation in a 75% full cylinder after (a) 
2 s, (b) 180 s and (c) 960 s of rotations. Color codes: brighter colors indicate higher poppy concentration and darker colors 
indicate higher millet concentration. FOV x x y: 6 cm x 6 cm. 



abrupt drop of the signal intensity near the cylinder wall, 
i.e., r = 30 mm. Further rotations lead to the formation 
of an annulus of poppy seeds, which is visible at 180 s, 
trapping some of the millet seeds at the core center, which 
corresponds to approximately 50% poppy by volume, in 
Fig. [6](b) and a trough in the ID profile in Fig. [71 at the 
core center. The maxima in all ID profiles correspond to 
the most poppy-rich region inside the cylinder. 

From this point in time, the behaviors of the 75% and 
82% full cylinders diverge. In the case of the 75% full 
cylinder, the maximum of the ID profile is only tempo- 
rary and disappears as radial segregation progresses to 
its end. This can be clearly seen in Fig.[6](c) in which the 
poppy and millet mixture in the core in Fig.EJb) has dis- 
appeared before the 960 s has elapsed. Equivalently, Fig. 
[7Ja) shows that this poppy-rich region initially appears 
at 8 mm away from the core center after a few seconds of 
rotations but gradually migrates towards the core center 
until they meet after 250 s of rotations. Subsequently, 
poppy seeds at the core center keep accumulating until 



its concentration reaches 100%. In the case of the 82% 
full system, however, the maximum of the ID profile is 
clearly permanent as it persists through to the comple- 
tion of radial segregation, which is also the onset of axial 
segregation after 2100 s of rotation as shown in Fig. [71(b). 

This formation and disappearance of the core center 
during radial segregation are not only interesting dy- 
namic behaviors in themselves but also tell us some in- 
formation about the avalanche layer thickness. 

Firstly, for a 75% full cylinder, the disappearance of 
the core center means that the avalanche layer extends 
well into the center of rotation. From geometrical con- 
siderations, this distance (i.e., from the top of the free 
surface to the center of rotation) is 12 mm, which means 
that the avalanche layer in a 75% full cylinder, rotated 
at 35 rpm is at least 12 mm. In the literature to date, 
there is a range of reported values for the free surface 
layer thickness. One author reported a thickness of 6 
mm for a pile of mustard seeds in an avalanche [22[, in 
which case the core center in our system is buried well 
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FIG. 7: ID concentration profile across the cylinder radius in (a) 75 and (b) 82% full cylinder during radial segregation at 35 
rpm. Zero distance from core center indicates the rotation axis whereas maximum distance means the cylinder wall. 



under the free surface. However, it has also been argued 
that the free surface layer may be thicker than previ- 
ously thought where depths of 13 mm at 22 rpm [23[ and 
18 mm at 30 rpm [241 ] . both in a cylinder of 70 mm in 
diameter, have been measured by MRI. 

Secondly, for an 82% full cylinder, the core is clearly 
permanent during radial segregation as mentioned above. 
The location of the final, stable maximum can be re- 
garded as the boundary between the diffusion-dominant 
core and the convection-dominant free surface region — 
8.5 mm from the core center. This corresponds to a dis- 
tance of 7.5 mm from the free surface to the boundary, 
which can be considered as the avalanche layer thickness. 
However, it has just been shown in the 75% full system 
that the distance from the core center to the free sur- 
face is 12 mm. This could be attributed to two reasons. 
First, the avalanche layer thickness could be dependent 
on the filling fraction. In the case of an 82% full, the 
free surface area is much smaller than that of a 72%, so 
that the distance available for acceleration of seeds from 
the top to the middle of the avalanche layer is smaller, 
thus the avalanche layer is dimmed thinner. Second, the 
avalanche layer could have significant effect on the core 
underneath so that in a 75% cylinder, the millet seeds ini- 
tially trapped in the core center have diffused out, leaving 
a more or less pure poppy core. 

As for the role of core diffusion, we examine the poppy 
concentration at the core center during radial segrega- 
tion. In order to separate the effect of the avalanche 
layer, we only look at the core center in the 82% full 
cylinder, which has been suggested above to be buried 
well under the avalanche layer, i.e. 8.5 mm. As shown in 
Fig. 0(b) , this concentration increases in a monotonous 
fashion from 25% by volume, i.e., homogeneously mixed 
state, to 67%. This indicates an outward flux of the mil- 
let seeds, which is most likely to be brought about by 
core diffusion as no convection is observed in this area. 



This confirms the importance of core diffusion, along side 
with the free surface, by which segregation takes place. 

To assess the relative significance of core diffusion and 
the free surface, we look into the accumulation rate of 
poppy seeds at the core center and the most poppy-rich 
region, i.e., at the maximum of the ID profile, as plot- 
ted in Fig. [8l As can be seen in Fig. [81(a) for the 75% 
full cylinder, the maximum and core concentrations are 
essentially the same except during the first 200 s when 
there still exists a millet and poppy mixture core. The 
poppy concentration in the core increases at an exponen- 
tially decaying rate as shown in the best fit equation. 
The rate at which the poppy seeds accumulate inside the 
core is represented by the exponential term of 0.01 s -1 . 
The 82% full system, on the other hand, has distinctive 
signal intensities at the core center and the most poppy- 
rich region, also defined as the boundary between the 
core and the avalanche layer, as shown in Fig. [51(b), con- 
firming the existence of the millet and poppy mixture in 
the core. It clearly proves the diffusion of millet seeds, 
in other words the accumulation of poppy seeds, inside 
the core at a rate represented by the exponential term of 
0.002 s _1 . In comparison with that of the 75% full sys- 
tem, i.e., 0.01 s _1 , it, however, occurs much more slowly 
due to the fact that the diffusion-dominant core is thicker 
for the 82% full system from geometrical considerations. 

The most interesting feature to note is that the ac- 
cumulative concentration at the poppy annulus follows 
a double exponential equation. Fitting the data to this 
double exponential equation yields a perfect fit whereas 
that to a single exponential equation is not satisfactory 
as shown in Fig. Etc). The fitted single exponential line 
shoots up very fast and levels out abruptly at t = 70 s 
while the experimental data levels out much more grad- 
ually. In the double exponential equation, the first term 
has the same exponential power of 0.01 s _1 as the 75% 
full system, but the second exponential power is higher, 
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FIG. 8: Poppy concentration at the core center and the maximum of the ID concentration profile across the cylinder radius in 
Fig. [7] during radial segregation in (a) 75% and (b) 82% full cylinder, (c) Comparison of the single and double exponential fit 
to the data for the maximum concentration in the 82% full cylinder. 



i.e., 0.14 s _1 , indicating a higher rate of poppy seed accu- 
mulation. This suggests two processes going simultane- 
ously at this boundary: diffusion inside the annulus and 
convection outside it from the avalanche layer. This is 
another confirmation that radial segregation is brought 
about by both convection in the free surface and diffusion 
inside the core. When the system is less than 75% full, 
all seeds inside the cylinder reach the avalanche layer in 
every rotation of the cylinder and thus models and the- 
ories based on the dynamics of the avalanche layer are 
appropriate. However, when the cylinder is full like those 
presented in this work, diffusion plays a significant role 
in the physics inside the core along the entire rotation 
axis. 

It should be again noted that all 3D structure images 
presented so far are acquired when the system is at rest. 
Radial segregation in an 82% full cylinder, rotated at 10 
rpm is also captured in real time by 2D FLASH MRI 
and the movie can be viewed at the provided link [25[. 
Some pictures extracted from this movie, showing initial 
radial segregation within the first 13.8 s of rotation, can 
be seen in Fig. [9j It is evident that initial formation 
of the poppy core during radial segregation within a 
few rotations is due to the stopping and percolation of 
poppy seeds, as suggested by Ref. [8[. 

Radial segregation: 82% full cylinder at 10 rpm: [25[ 



B. Axial band formation 

We present 3D MRI data for sections of a 50% full, 50 
cm long drum, rotated at 35 rpm. The progress of the 
band formation is captured, thus its rate can be calcu- 
lated. Axial band formation is a familiar topic but we 
demonstrate in this section that full 3D structure of the 
cylinder can be captured by MRI, thus is useful for mod- 
eling work. We will also address the role of core diffusion. 
For this purpose, we will present MRI data for an 82% 
full system. 



Figure [10] shows the progress of axial segregation in a 
22 cm section near one end of the rotating drum, which 
is made by putting together two FOV 12 cm images end- 
to-end (overlapping 2 cm). There is a small band right 
at the end cap (toward the right of the images at y = 
0), which is not shown but we capture the formation of 
two real bands adjacent to it. At 270 s, the bands start 
to form with poppy seeds breaking off at two locations 
along the rotation axis. At 899 s, the formation of two 
axial bands is more or less complete and it is evident that 
there are poppy seeds along the entire axis, ranging from 
40% to 60% by volume as indicated by the red and green 
colors, connecting the two poppy bands. 

Similar to the procedure employed for analyzing radial 
segregation, a series of 3D images like those in Fig. [10] 
are acquired and collapsed into ID profile showing the 
concentration of poppy seeds along the rotation axis by 
averaging the signal intensity along the x and z axes due 
to their symmetry. Results are shown in Fig. [TTJ Dur- 
ing the course of band formation, the band position does 
not change, but only shortens in width and increases in 
height, i.e., poppy concentration. Axial band formation 
therefore involves essentially the diffusion of millet seeds 
from the nascent poppy band to the adjacent nascent mil- 
let band, which is in agreement with the mechanis m by 
the dynamic angle of repose as proposed by Kakalios [2l| . 
The rate at which these two bands are formed can be de- 
termined by fitting an equation to the bands' maximum 
poppy concentration. An increase of poppy accumula- 
tion at an exponentially decaying rate similar to Fig. [8] 
is observed. 

It is of interest to examine if the dynamics of the 
avalanche layer is the only mechanism of axial segrega- 
tion, as claimed in many studies. For example, an at- 
tempt to verify this mechanism was made by conducting 
a simple experiment, in which a co- axial cylinder was in- 
serted into the segregating cylinder and axial segregation 
was no longer observed [26[. The author then concluded 
that free surface was the only mechanism for segrega- 
tion since the experiment eliminated free motion along 
the surface. However, this setup also removed the oppor- 
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FIG. 9: 2D FLASH pictures capturing radial segregation in an 82% full cylinder, rotated at 10 rpm during the first 12 s of 
rotation. The timings of these snap shots from left to right, top to bottom are: 0.9, 1.8, 3.0, 4.2, 5.4, 6.6, 9.0 and 13.8 s. Color 
codes: blue: millet, red: mixture, yellow-white: poppy. 




FIG. 10: 3D structures showing the formation of two axial bands in a section near one end of a 50% full cylinder after (a) 270 
s and (b) 899 s of rotations at 35 rpm. 



tunity for core formation; hence a conclusion cannot be 
drawn about the sole mechanism of segregation. Instead, 
we conduct a different experiment, in which a rather full 
cylinder is employed. If the free surface is the only mech- 
anism, it is anticipated that axial segregation should not 
occur in such a full system since there is only a small 
amount of free surface and more importantly, there exists 
a core which is buried under the free surface as presented 
previously when discussing radial segregation. The cylin- 
der to be studied is the same as that for the radial segre- 
gation, which is 12 cm long and filled to 82% by volume. 
The final state of this system after 55,000 revolutions (26 
hours) is shown in Fig. [T2J The system still segregates 
like a 50% system so that there is one poppy band next to 
the left end cap, and the rest is filled with millet seeds. 
As we saw earlier, there exists a core along the entire 



length of the rotation axis with the diameter of at least 
7.5 mm buried under the avalanche layer. However, this 
core only remains inside the poppy band and has com- 
pletely disappeared from inside the millet band. This is 
a clear indication of core diffusion leading to complete 
axial segregation, which is in line with the sub-surface 
axial segregation observed previously [12]. This core of 
millet and poppy mixture is not visible to the eye or by 
ordinary measurements but can be detected by MRI. As 
shown in the cut-through view in Fig. [12] there is a core 
of ~40% poppy, located on the rotation axis but inside 
the poppy band. 

Real time movies showing the band formation can also 
be acquired by 2D FLASH MRI. A 50% full and 12 cm 
long system is chosen so as to keep the segregation time 
sufficiently short for data acquisition. The movie J2?J 
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FIG. 11: Axial band formation in a 22 cm section of a 50% 
full and 50 cm long cylinder, rotated at 35 rpm. 




FIG. 12: 3D structure of the final state of axial segregation 
in an 82% full and 11 cm long cylinder, rotated at 35 rpm. 

is recorded over 4.7 min of axial segregation and some 
intermediate pictures are shown in Fig.[13l In this movie, 
one can clearly see the formation of two poppy bands 
at the two end caps and one band in the middle of the 
cylinder. This is just to demonstrate the applicability 
of MRI for studying the dynamics of granular materi- 
als and that it can provide useful data for modeling work. 

Axial band formation: 12 cm long and 50% full cylin- 
der: 



C. Axial band traveling and merging 

Upon "full" axial segregation, bands have been ob- 
served to travel and merge over a long time scale (many 
hours) by some authors [H, [13, EE O- This phe- 




FIG. 13: Axial band formation in a 50% full and 12 cm long 
cylinder, rotated at 35 rpm during the first 2.5 minutes of 
rotation. The timings of these snap shots from the top to 
bottom are: 0.3, 3.6, 18.8, 95.0, 133.7 and 284.4 s (4.7 min). 
Color codes: blue: millet, red: mixture, yellow- white: poppy. 

nomenon is observed in our study on both short (within 
13.5 min) and long (after 33 hours) time scales in a series 
of successive merging events. In order to keep the time 
scale sufficiently short for data acquisition, only 50% full 
systems are considered. The long time scale band trav- 
eling and merging are first studied in a system of a 50 
cm long rotating drum by 3D stationary RARE. Such a 
long cylinder is chosen so as to permit the formation of 
many axial bands (8 poppy bands in this case), followed 
by band traveling and merging events. The real-time 2D 
FLASH movies of both band traveling and band merg- 
ing in a 12 cm long cylinder are subsequently presented. 
This short cylinder length is chosen for 2D FLASH so 
that both band traveling and merging occur within 13.5 
min of rotation and are feasible to produce a movie in- 
stead of over 33 hours. 

As mentioned previously, there are 8 bands of poppy 
seeds formed at equal spacings of 7 cm at full axial seg- 
regation in a 50% full and 50 cm long drum. Further 
rotations over the next 70,000 revolutions or 33 hours 
result in these bands to merge pair-wise until there are 
only two bands in the entire cylinder, one is the small 
original band at one end of the cylinder and the other 
is a long band located at the other end. The number 
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FIG. 14: In a 50% full and 50 cm long cylinder: (a) number of 
bands over time due to band merging and (b) band traveling 
prior to two separate merging events. 



of poppy bands over time follows a logarithmic decaying 
rate, as shown in Fig. Ha), which was also observed by 
Ref. 0. 

There are many merging events, which are always pre- 
ceded by band traveling, starting from 8 poppy bands 
at full axial segregation, only two of which are chosen 
as the case study for band traveling. It is observed that 
prior to each merging event, mainly one band is "ac- 
tive" and travels gradually toward the other band, which 
remains essentially stationary. During traveling, bands 
maintain their shape (band with side bulges as shown in 
Fig.fTOlfb)) and size. FigurefTWb) shows the displacement 
of the traveling bands in the two separate merging events 
mentioned above over time. The linearity of both curves 
in Fig. [T4lf b) indicates a reasonably constant, as previ- 
ously reported [lj] , but very sluggish traveling speed of 3 
± 0.2 /im s _1 (taken from the gradient of the two curves). 
However, once the bands are adequately close — approx- 
imately 3 cm apart — band merging occurs much more 
rapidly. 

Unlike band traveling, poppy bands remain stationary 
during merging but change their shape by the movement 
of millet seeds through the poppy bands. The merging 
of poppy bands can be seen in Fig. [15] and summarized 
as follows. Poppy seeds first travel along the rotation 
axis from the existing bands, forming a thin line along 
it, connecting the two bands. This can be seen in Fig. 
[T5lfa) where the area along the rotation axis between the 
bands appears green, indicating a poppy and millet mix- 
ture consisting of 30-60% poppy by volume. Millet seeds 
then keep traveling out of this area through the poppy 
bands over many stages such as one in Fig. fTST b) until 
the concentration of poppy seeds in the bands and the 
area between are equal as shown in Fig. [T5lc). A series 
of images like those are recorded to confirm the merging 
mechanism. 

These images confirm a similar claim in Ref. [l2| that 
the merging of bands is not simply their traveling until 
they meet each other, but involves the removal of mil- 
let seeds from between the poppy bands. In that paper, 
the camera technique was used so only observations on 
the surface could be made. However, in this paper, full 
3D images of particle distribution are available, show- 
ing the same process happening in the bulk materials. 



Also, useful information listed below, such as the merg- 
ing mechanism and rate, can be extracted: 

First, the images confirm the diffusion of millet seeds 
through the poppy bands in that the poppy concentration 
in the poppy bands is slightly reduced during merging as 
shown in Fig. [T6la). This indicates the presence of lower 
intensity millet seeds inside the poppy bands during the 
merging process. This is also in agreement with visual 
inspection of the presence of millets, although mostly on 
the surface of the poppy bands. 

Second, it is evident that band merging involves the 
transport of millet seeds from between the poppy bands. 
Figure [T6la) shows the concentration average along the 
cylinder length of the last merging event, in which the 
band on the right hand side is an active band and the 
other is a stationary band at one end of the cylinder. Af- 
ter 11,700 s (3.25 hours) of rotation, these bands start to 
merge. The trough between the bands is initially pure 
millet seeds, i.e., 0% poppy. The subsequent movement 
of millet seeds out of the area between the two merging 
poppy bands is illustrated by the increase of poppy con- 
centration in this region until it equates that in the poppy 
bands themselves. At this stage, there is still a consid- 
erable amount of millets seeds — 37% by volume. Mil- 
let seeds then keep moving outward, leaving the merged 
poppy band more concentrated until it reaches the con- 
centration of a fully developed band. 

Third, the rate at which band merging occurs can be 
derived from the acquired images, as shown in Fig. [T6lf b). 
which shows the concentration minima in the area be- 
tween the two merging poppy bands over time. It is 
worth noting that this graph primarily focuses on the 
time interval in which fast merging occurs. Therefore, it 
does not show sluggish processes of band traveling prior 
to merging not the millet removal process in which the 
poppy concentration gradually increases after merging. 
A logistic trend is recognized, reinforcing our claim that 
band traveling occurs very gradually but once merging 
starts, it takes place very rapidly. Similar graphs, though 
not shown here, can be generated with two other merging 
events prior to the one examined above. More abrupt 
logistics curves are obtained, which means merging oc- 
curs more rapidly. This could be attributed to the fact 
that the second and third last merging events involve two 
active bands, meaning bands away from the end cap. 
Therefore millet seeds are free to move through both 
poppy bands during the merging process hence speeds 
up merging. However, the last merge is between a fixed 
band at the end cap and an active poppy band so that 
millet seeds in between the bands have only one option of 
traveling through the one active band, thus slows down 
the millet transport process. 

It should also be noted that the "final" poppy bands 
are in fact not pure poppy seeds but contain millet seeds. 
By MRI, we can quantify this accurately to be 60% poppy 
by volume. As mentioned, the newly merged poppy band 
still contains considerably amount of millet seeds, i.e., 
37% poppy by volume, but becomes more concentrated 
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FIG. 15: 3D structures showing the progress of one band merging event in a section of a 50% full and 50 cm long cylinder at 
(a) 67790 s, (b) 69711 s and (c) 72301 s. The left hand side band is a wall band and therefore stationary during merging. 




FIG. 16: (a) ID poppy concentration average along the rotation axis of the band merging in Fig. \TE\ over time, (b) Poppy 
concentration of the trough in (a) over time. 
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subject to further rotations until it reaches 60%, though 
not shown in Fig. [T6lfb). A confirmation of this figure 
is by visual inspection and using the length of the final 
poppy band to compare with the known volume of seeds 
in the system. The same poppy concentration is found. 
Also, the same concentration is calculated for all inter- 
mediate poppy bands over the course of 80,000 revolu- 
tions (38 hours). This is, therefore, considered to be the 
ultimate poppy concentration. To our knowledge, only 
visual inspection has been reported that there are some 
millet seeds in the poppy bands but no quantification is 
given due to the lack of appropriate tools. This is the 
first time the poppy concentration inside poppy bands is 
quantified accurately and non-intrusively by MRI. 

As mentioned previously, band traveling and band 
merging are also observed on a much shorter time 
scale. A cylinder of 12 cm in length is chosen, in which 
three temporary poppy bands are formed at full axial 
segregation. This length is short enough to be fitted 
within the FOV of the magnet but also long enough to 
allow three poppy bands to form and merge eventually. 
The middle band then travels toward and merges with 
one band at an end cap. For the first time, this dynamic 
process is imaged in real time by 2D FLASH movie as 
shown at the provided link and Fig. [T71 capturing 13.5 
min of rotation. 

Band traveling and merging: 12 cm long and 50% full 
cylinder: [H[ 

IV. CONCLUSIONS 

In this paper, we present full 3D MRI data of radial 
and axial segregation in rotating drums. 2D FLASH real 
time movies are also recorded showing the mechanism of 
initial radial segregation, the formation of axial bands as 
well as band traveling and merging. These data should be 
very useful for modeling work as all spatial information 
at all intermediate stages of radial and axial segregation 
are captured by MRI. 

Some new observations are made in this work. First, 
during radial segregation, there is an inner core formed 
along the rotation axis, which eventually disappears in a 
75% full cylinder but is permanent in an 82% full cylin- 



der. From this, the avalanche layer thickness can be in- 
ferred. Second, band traveling and merging are observed. 
Prior to each merging event, mostly one poppy band trav- 
els at an almost constant speed of 3 fim s _1 , until the two 
merging bands are approximately 3 cm apart. The merg- 
ing process then involves the sideway traveling in both 
directions of millet seeds from between the poppy bands 
until the concentration here equals that of the two poppy 
bands themselves. Millet seeds then keep traveling out of 
the merged poppy band, thus shortening this band, un- 
til the ultimate poppy concentration of 60% by volume 
is reached. The speed of merging follows a logistic rate, 
which is more rapid if both poppy bands are active, i.e., 
away from the end cap. 

Initial radial segregation (within the first few rota- 
tions) is due to the stopping and percolation of the poppy 
seeds. Subsequently, the mechanisms of segregation are 
identified to be both free surface and core diffusion. The 
latter has largely been ignored in the literature but is 
shown in this work to become significant in cylinders of 
more than 75% full since there exists an inner core buried 
under the avalanche layer. This core can be clearly seen 
during radial segregation but becomes axially segregated 
during axial band formation in an 82% full cylinder. A 
future work will follow the current one to study the sub- 
diffusive nature of this core diffusion phenomenon. 

In addition, it is demonstrated that MRI is a quanti- 
tative technique in which the seed concentration at all 
spatial points can be directly calculated from the sig- 
nal intensity. The concentration of all poppy bands, in- 
cluding all initial, intermediate and final poppy bands, 
is quantified accurately for the first time to be 60% by 
volume of poppy seeds. 
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